DIRAC ELECTRON IN THE ELECTRIC DIPOLE 

FIELD 

V.I.MATVEEV, M.M.MUSAKHANOV 
Department of Theoretical Physics, Tashkent State University 
700095, Vuzgorodok, Tashkent, Uzbekistan 
e-mail : yousuf @ aphi . silk. glas . ape . org 
and 

D.U.MATRASULOV 
Thermal Physics Department of Uzbek Academy of Sciences 
700135, ul.Katartal 28, Tashkent, Uzbekistan 

February 1, 2008 

Abstract 

Dirac equation for the finite dipole potential is solved by the method of the join of 
the asymptotics. The formulas for the near continuum state energy term of a relativistic 
electron-dipole system are obtained analytically. Two cases are considered: Z < 137 and 
Z > 137 

Dirac equation for Z > 137 is solved by usual method of cut-off potential at small 
distances. 







We would like to discuss the problem of the Dirac equation for an electron in an electric 
dipole field which is closely related to the problem of exotic heavy quarkonium states too. We 
mean the problem of the motion of light valence and vacuum quarks in the color field of heavy 
quark- ant iquark pair. The motion of an electron under the influence of an electric dipole one of 
the classical problem of nonrelativistic quantum mechanics. Energy eigenvalues of such system 
was found by many authors see e.g. ]T|, |2], Several authors have calculated the critical value 
of the dipole moment at which energy term reaches the boundary of continuum [l], 0. 

Solution of corresponding relativistic problem is more difficult mathematical problem, since 
variables of the Dirac equation can't be separated at any orthogonal system of coordinates. 
Near continuum energy state of the relativistic electron, moving in the field of finite electric 
dipole, we shall find by the method of the join of the asymptotics which successfully was 
applied to solution of two-center Coulomb problem for the Dirac equation (with same signs of 
charges) || [7|. Previously this problem was considered in our paper |§. 

In the present paper we also use asymptotical separability of variables at large distances in 
the relativistic two-center Coulomb problem and additionally discuss the nonrelativistic limit 
of our results. Then the squared Dirac equation takes a form of the Schrodinger equation at 
which variables can be separated. Let us consider relativistic electron moving in the field of a 
finite electric dipole which is composed of charges +Za and — Za are separated by distance R. 
The potential of this system is given by 

y= Za_Za 
n r 2 

where r, is the distance from i th charge, a -1 = 137 (the system of units h — m — c — 1 
is used here and below). The motion of this electron is described by Dirac equation which can 
be written as 

< 2 > 

where H = ap + (3 + V is the Dirac Hamiltonian, a and f3 are the Dirac matrices. 

In the spinor representation stationar squared equation can be written in the components 

(i'A 

ipi and ip2 of a byspinor ip — \ as in [5] 

[(e- Vf + A + iffW-l]^i (3) 

[(e - V) 2 + A + iaVV - 1]^ 2 (4) 
where e is energy of an electron, a are standard Pauli matrices. 
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At large distances from dipole we can neglect by terms V 2 and ia\/V with respect to V 
hence for each component we have equation 

[A - 2eV + e 2 - l]ip = (5) 

Eq.(||) is the Schrodinger equation with potential eV and with energy — \ (e 2 - 1) So, at 
large distances the wave function has separated form 

^ = (£2?i)i/ 2 (i Z^yf2 exp(im>p) (6) 

where £, rj, ip are the spheroidal coordinates which are defined by £ = (t*i + r2)R~ 1 , r\ = 
(ji — r 2 )R~ 1 , f = arctg(^), where y and x are the cartesian coordinates, R is the distance 
between charges, 1 < £ < oo , — 1 < 77 < 1 , < ip < 2n 
The radial and angular equations become 

u"(o + w + ^ + ^r^f ]m = (7) 



v "(v) + l-p 2 + ^—r + Zv Wv) 



Dr] — A 1 — m 2 
f] 2 (1 — I] 2 

where p 2 = — ^-(s 2 — 1), D = 2 | e \ RZa, A is the constant of separation. 
The boundary conditions are 

U{1) = 0, oo — > 0, 7(±1) = (9) 

Asymptotes of eq.(0) for large £ is [1] U(£) = (pQ^KiuipOi where K iu is the McDonald 
function, v 2 = A — |. Solution of eq.(§) in the vicinity of 77 = was also found in ref. [1]. 
This asymptotes has the form 

V(rj) = [z'(r])]^ Ai(-Dh) (10) 

where Ai is the Airey function, z(rj) = [| $ ^^(l — r? 2 ) -1 / 2 ^]! and 2/(77) = 

Thus for large p£ and for small 77 asymptotes of wave function can be written in the form(for 

771=0) 

ip = [ — ^ J ^(K)% 2 - lj'^in^/nf— + arfl7(l + ii/)]) (11) 
Here we have used asymptotes of function K iv {x) for x — > [|J 

K iu (x) — ► [ — 1 2 (p£) 2sini/ln[ — h arg / -f(l + ii/)] 

ushTru x 

and the asymptotes of Airey function [|15[ Ai(z) — > C\ — c 2 z for z — > 0, where c\ = 0.355, c 2 = 
0.259, and the relation 2/(77) « 1 for 77 — > 0. 



Following by Popov || as an asymptotes of wave function at small distances we will take the 
product of relativistic one-center wave functions. Here we must consider two cases: Z < 137 
and Z > 137. 

Let Z < 137. In this case asymptotics (at small distances from center) is given by [5] 

0! = r~ 1+7 (12) 



where 7 = (1 — Z a" 

Thus asymptotes of two center wave function at small distances from dipole can be written 
as follows: 

$x~{Tir*y<- 1 = {e-7?y<- 1 (13) 

Let now Z > 137. In this case energy becomes imaginary ("fall to the center") || [T0[| . In 
order to solve the Dirac equation in this case we must cut-off the Coulomb potential at small 



distances lO, |12 



General form of the cut-off potential is |I~3| , |TJ| 

V{r) = < 



Za 
r 1 



for r > b 
¥/(r) , for < r < b 



where b is the radius of the nuclei. 



For f(r) = l (surface distribution of charge) asymptotes of wave function has the form||T3" 



02 ~ r (14) 

Hence for the case of Z > 137 asymptotics of two center wave function at small distances from 
dipole can be written in the form 

ip ~ rir 2 ~ £ 2 — rj 2 . (15) 

So in general case one can write for the asymptotics of wave function at small distances from 
dipole (for R « 1) 

~ (e - v 2 f, (i6) 

where 

( 7 - 1 , for Z > 137 



For small 77 



1 , for Z> 137 

fa ~ i 2P (17) 
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Since | e |~ 1 asymptotes of wave function at small and at large distances join with other in the 
region 1 << £ << K Equating logarithmic derivatives of functions (pj]) and (|TT|) and taking 
into account that in this region ln-^ ~ In- we have 

1 2 

- + vctg[vln- + argT(l + iv)) = -2(3 (18) 
Z p 

From this expression we get 

p = exp(—argT(l + iv) arcctgl ———] — 2ln2) (19) 

v v 2v 

In order to find the e(R) from (|T^) we must know the parametr v. Here we use results of 
ref.0 where dipole moment was found as a function of v. In application to our problem for 
ground state this formula can be written as 

D=-^i[l + ^ (V + 6^ 2 -l)] (20) 

where D = 2 | e \ ZRa. From this equation we find v : 

v = [^[D - 2 | e | ZRa - D cr ] - ^-(e 2 - l)]h (21) 

where D cr = ^[1 - £]. 

is the critical value of dipole moment at which term e(R) reaches the boundary of continuum. 
Taking into account relation between p and e we find 

e ± = ±[l-^exp(u(v,Z))} (22) 



where 



2 2 1 + 4/? 

w(v, Z) = —argTil + iv) arcctgl ] — 2/n2 (23) 

v v 2v 



Thus for Z < 137 we have transcendental equation for e(R) 

2 r 2 n 2 3-47 , . 

e± = ±1 =p —exp[—argT(l + iv) arcctg — 2/n2 (24) 

R 2 v v 2v 

where v defined by (21). 

For Z < 137\/7/4 eq.(p4]) can be solved by iterations (since iteration procedure converges 
only for 3 — 47 < ) So for the first iteration we have 

+i^ 2 r r2 (i) , r 2 (i)(3-4 7 ) 

e± = ±1 =F ^-exp T-arcctq r (25) 

R 2 L 7i{AZaR-2D cr ) 2 4n{4ZaR - 2D cr ) * 

Formula (p5|) is the near continuum state energy term of an electron in the electric dipole 
field. As well known nonrelativistic limit gets near the continuum. Since for near continuum 



states AZaR — 2D cr ~ 0, arcctg — > 7r. Thus in order to obtain nonrelativistic limit we must 
do following replacements: 1)7 — > 1, 2) arcctg — > result we have 

e± = ±1 =F T^h U ] (26) 

71(4^/? - 2D cr .)2 

This formula for e± coincides with known formula from || for the nonrelativistic electron-dipole 
system. By expanding arccotangent into series in formula ( ^5|) we calculate the correction to 
eq.@ 

2 r 2 fi 4 

£ ~ T^expl ^ r - — ] (27) 

R 2 1 n(4ZaR-2D cr ) l 2 3 - 4 7 J v ; 

From this formula it is easy to see that nonrelativistic limit take place for (AZaR — 2D cr )^ « 
0.8. We note that this nonrelativistic limit can be also obtained directly from eq. (|i~9"|) . For 
Z > 137 we have equation 

2 2 2 5 2n 

e± = ±1 =p — ^exp[—argT(l + iv) arcctg 2ln2 ] (28) 

R 2 v v 2v v 

Solving this equation by iterations for the first iteration we have 

2 r 2 fi) 5r 2 (i) 

£± ~ ±1 =F ~^exp[ ^ T {arcctg ^ j )] (29) 

R 2 1 7c{AZaR-2D cr ) L 2 K 47r{4ZaR - 2D cr ) l 2 - n 

Derived analytical formulas have to be useful for further numerical calculations in nonasymp- 
totical region. 

We are planning to generalize our results for the Coulomb plus confinement potential case 
for the calculations of the properties of the exotic quarkonium states. 
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